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1.  INTRODUCTION 

A  critical  step  in  the  practice  of  ground-based  optical  inter¬ 
ferometry  is  the  measurement  of  suitable  calibrators  (i.e.,  Mo- 
zurkewich  et  al.  2003).  These  allow  compensation  for  various 
instrumental  and  atmospheric  effects  that  reduce  the  measured 
fringe  contrast.  The  ideal  calibrator  is  bright  (short  of  inducing 
nonlinearities  in  the  detector)  and  has  either  an  accurately 
known  visibility  or  a  diameter  that  can  be  reliably  estimated 
to  be  small.  In  the  latter  case,  stars  that  are  bright  and  have 
small  angular  diameters  are  necessarily  high-surface-brightness 
objects  (i.e.,  early-type  stars). 

Early-type  stars  bring  their  own  problems.  In  practice,  the 
need  to  have  calibrators  relatively  close  to  the  science  targets 
results  in  objects  being  chosen  that  have  diameters  that  are 
larger  than  optimal.  Unfortunately,  as  has  been  known  almost 
since  the  inception  of  spectroscopy,  stars  earlier  than  F5  tend 
to  rotate  rapidly.  This  leads  to  oblate  figures  and  significant 
gradients  in  surface  temperature  (von  Zeipel  1924),  resulting 
in  visibilities  that  deviate  from  predictions  by  more  than  the 
error  estimates,  and  even  worse,  visibilities  that  can  change 
significantly  with  hour  angle.  Reflecting  this  fact,  many  lists 
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of  calibrators  (e.g.,  Borde  et  al.  2002;  Merand  et  al.  2005)  have 
excluded  hot  stars  for  the  most  part. 

The  actual  rotation  state  of  individual  stars  is  rarely  known, 
since  the  spectra  give  only  the  projected  equatorial  velocities 
V  sin  i,  usually  interpreted  only  through  line  widths.  Neverthe¬ 
less,  projected  velocities  can  provide  useful  constraints  on  the 
range  of  the  effects  on  the  visibilities.  In  this  paper,  we  use 
the  theory  developed  by  von  Zeipel  (1924)  to  estimate  the 
visibilities  over  the  range  of  parameters  consistent  with  a  star’ s 
observed  magnitude,  color,  projected  rotational  velocity,  and 
parallax,  and  its  location  in  appropriate  evolutionary  grids. 

In  the  next  section,  we  describe  the  models  and  discuss  their 
accuracy  by  comparing  the  diameters  deduced  from  our  cal¬ 
culated  visibilities  with  measurements  found  in  the  literature. 
Furthermore,  we  have  found  that  in  many  cases,  the  existence 
of  a  lower  limit  to  the  (polar)  surface  gravity  imposed  by  the 
zero-age  main  sequence  (ZAMS)  provides  an  important  con¬ 
straint  on  the  range  of  possible  inclinations,  and  in  §  3  we 
describe  its  impact  on  the  simulations.  In  §  4,  we  discuss  the 
results  for  some  representative  examples. 

We  have  compiled  lists  of  early-type  stars  that  might  be  used 
as  calibrators.  Using  typical  parameters  for  some  currently  op¬ 
erating  interferometers,  we  characterize  for  each  object  the 
range  of  visibilities  that  might  be  induced  by  rotation.  In  §  5, 
we  describe  the  selection  of  these  objects  and  the  format  of 
the  catalogs. 

A  preliminary  report  of  this  research  has  recently  appeared 
(Yoon  et  al.  2006),  in  which  we  focused  on  the  more  practical 
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aspects  of  rotation  and  the  calibration  problem.  Some  details 
that  we  skip  here  can  be  found  in  that  report. 

2.  THE  EARLY-TYPE  STARS  AS  CALIBRATORS 

2.1.  The  Theory  of  Rotating  Stars 

The  theory  of  rotating  stars  was  studied  by  von  Zeipel  in 
1924.  He  showed  that  stars  in  solid-body  rotation  will  adopt 
the  figure  of  a  Roche  spheroid  to  first  order.  He  further  showed 
that  if  radiative  equilibrium  holds,  the  emergent  flux  varies 
over  the  surface  proportionally  to  effective  gravity;  that  is, 
effective  temperature  is  proportional  to  the  fourth  root  of  the 
local  effective  gravity.  This  theory  has  been  tested  not  only  in 
close  binary  systems,  but  also  through  interferometric  obser¬ 
vations  of  isolated  stars  (e.g.,  Domiciano  de  Souza  et  al.  2005; 
Peterson  et  al.  2006a,  2006b;  van  Belle  et  al.  2006). 

2.2.  The  Visibility  Calculations 

The  calculation  of  the  visibilities  was  carried  out  as  described 
in  Peterson  et  al.  (2006a).  Briefly,  given  values  for  the  Roche 
parameters  (e.g.,  §  3.1),  the  position  angle  of  the  rotational 
axis,  and  the  orientation  of  the  projected  interferometer  base¬ 
line,  we  construct  a  grid  of  points  across  the  projected  disk  of 
the  star,  with  axes  parallel  and  perpendicular  to  the  projected 
baseline. 

Following  the  developments  in  Collins  (1963)  and  Hardorp 
&  Strittmatter  (1968),  we  calculate  local  surface  values  of  ef¬ 
fective  temperature,  effective  gravity,  and  the  direction  of  the 
normal  at  the  grid  points,  and  using  the  Kurucz  (1993)  models 
as  parameterized  by  Van  Hamme  (1993),  the  specific  intensity 
along  the  viewing  angle.  Intensities  perpendicular  to  the  base¬ 
line  are  summed,  creating  the  “strip  brightness  distribution” 
(SBD).  The  Fourier  transform  of  each  pixel  along  the  baseline, 
taken  as  a  hat  function,  is  calculated.  The  result  is  weighted 
according  to  the  normalized  SBD,  multiplied  by  a  phasor  ac¬ 
cording  to  the  “shift  theorem”  (  Bracewell  1965),  and  summed 
to  produce  an  approximate  complex  visibility  V  (Born  &  Wolf 
1999).  The  squared  visibility  is  then  simply  =  \^.  We  do 
not  treat  the  related  issue  of  the  problems  created  for  phase 
measurements  due  to  rotation  of  the  calibrators  (Peterson  et  al. 
2006a). 

2.3.  The  Visibility  Calibration 

In  what  follows,  we  calculate  the  variations  of  the  visibilities 
that  are  expected  for  a  given  star  over  the  possible  range  of 
rotation  parameters.  The  primary  purpose  of  these  calculations 
is  to  characterize  the  range  of  visibilities,  but  it  is  also  of  interest 
to  know  how  reliable  the  calculated  visibilities  are  themselves. 

A  good  check  is  to  compare  the  angular  diameters  we  would 
predict  in  the  absence  of  rotation  with  the  measured  diameters 
reported  in  the  literature.  To  this  end,  we  list  in  Table  1  the 
long-baseline  angular  diameters  (0„eas)  for  all  stars  listed  in  the 
Catalog  of  High  Angular  Resolution  Measurements  (Richichi 


et  al.  2005)  in  this  color  range,  taking  care  to  eliminate  multiple 
references  to  the  same  measurements.  Since  our  calculated  di¬ 
ameters  (^caic)  implicitly  include  limb-darkening,  in  the  few 
cases  in  which  only  uniform  disk  diameters  are  quoted,  it  is 
necessary  to  correct  for  limb-darkening  effects.  This  was  done 
with  the  simple  device  given  by  Hanbury  Brown  et  al.  (1974), 
using  linear  limb-darkening  coefficients  appropriate  for  the 
star’s  temperature  and  the  bandpass  of  the  observations  from 
Van  Hamme  (1993). 

The  comparison  between  our  photometric-based  calculated 
diameters  and  the  measurements  is  shown  in  Figure  1.  Sys¬ 
tematic  differences  are  quite  small.  We  have  solved  separately 
for  possible  additive  (^^eas  =  ^caic  +  o)  or  multiplicative 
(^meas  =  ^  ^caic)  offsBts.  The  solid  line  in  Figure  1  shows  the 
best  least-squares  fit  for  the  additive  case,  a  =  0.084  ± 
0.026  mas.  The  result  for  the  multiplicative  offset,  b  = 
1.019  ±  0.007,  which  for  clarity  we  do  not  show,  is  compa¬ 
rably  small.  We  conclude  that  our  visibility  calculations  are 
systematically  reliable  for  stars  in  this  spectral  region. 

On  the  other  hand,  it  is  evident  that  the  scatter  is  fairly  large 
in  Figure  1,  and  indeed  we  deduce  reduced  x^’s  of  12.15  and 
13.93  for  the  additive  and  multiplicative  cases,  respectively. 
This  we  attribute  primarily  to  the  effects  of  rotation,  which  are 
explicitly  not  modeled  in  the  calculations  or  measurements. 
This  explanation  can  easily  accommodate  (3  Cas  and  a  Aql. 
The  large  residuals  for  a  CMi  (Procyon;  at  F5V,  the  coolest 
object  in  this  list)  are  a  puzzle  and  may  originate  in  the  effects 
of  shallow  convection,  as  discussed  by  Robinson  et  al.  (2005) 
and  Aufdenberg  et  al.  (2005). 

3.  SIMULATIONS  AND  INTERPRETATION 

In  order  to  consistently  estimate  the  visibilities  and  accu¬ 
rately  estimate  the  variations  induced  by  rotation,  we  cannot 
simply  use  angular  diameters  provided  by  the  usual  color-mag- 
nitude-diameter  relations  (i.e.,  Barnes  et  al.  1978,  hereafter 
BEM78;  Mozurkewich  et  al.  2003).  Instead,  we  must  calculate 
models  that  are  consistent  with  each  star’s  color  and  magnitude 
and  with  its  other  measured  properties  over  the  range  of  allowed 
inclinations  and  position  angles.  We  describe  those  calculations 
next. 

3.1.  The  Simulations 

Roche  models  require  several  parameters:  inclination  angle 
i,  mass  M,  angular  diameter  6  (or  equivalently,  the  radius  R 
through  the  parallax  p),  polar  effective  temperature  and 
the  dimensionless  angular  velocity  w  ( =  O/Og,  where  Og  is  the 
angular  velocity  at  breakup).  Starting  with  initial  guesses  of 
mass  and  temperature,  we  iterate  until  we  match  v  sin  i  by 
varying  the  angular  velocity  parameter  (w)  or  inclination  angle 
(i).  Within  this  iteration,  the  program  must  also  iterate  for  con¬ 
sistent  values  of  ^eft.g  and  Op  to  match  the  observed  V  mag¬ 
nitudes  and  B  —  V  colors.  With  values  for  w,  i,  Op,  par¬ 
allax,  and  an  assumed  mass,  we  can  calculate  the  polar  radius 
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TABLE  1 

Measured  Values  and  Calculated  Values  of  Angular  Diameter 


HR 

Name 

^calc 

f'm. 

as 

Bandpass 

References 

21  . 

11  &  Cas 

1.89 

2.12 

+ 

0.05 

740  nm 

1 

82 . 

27  p  And 

0.55 

0.66 

+ 

0.08* 

H,  K 

2 

269  . 

37  And 

0.67 

0.67 

+ 

0.11* 

H,  K 

2 

2421  . 

24  Y  Gem 

1.40 

1.39 

+ 

0.09 

443  nm 

3 

2491  . 

9  a  CMa 

6.12 

6.04 

+ 

0.02 

K 

4 

5.89 

+ 

0.16 

443  nm 

3 

5.99 

+ 

0.11 

800,  451  nm 

5 

2943  . 

10  a  CMi 

5.03 

5.45 

+ 

0.05 

K 

6 

5.45 

+ 

0.05 

451,  550,  800  nm 

5 

5.50 

+ 

0.17 

443  nm 

3 

6.44 

+ 

0.25 

706  nm 

7 

5.46 

+ 

0.08 

740  nm 

8 

4090  . 

30  LMi 

0.54 

0.59 

+ 

0.06* 

H,  K 

2 

4534  . 

94  ^  Leo 

1.34 

1.31 

+ 

0.09* 

443  nm 

3 

5447  . 

28  o  Boo 

0.69 

0.78 

+ 

0.04* 

H,  K 

9 

6556  . 

55  a  Oph 

1.54 

1.63 

+ 

0.13 

443  nm 

3 

6629  . 

62  y  Oph 

0.60 

0.65 

+ 

0.06* 

H,  K 

2 

7001  . 

3  a  Lyr 

3.20 

3.23 

+ 

0.03 

451,  500,  800  nm 

5 

3.24 

+ 

0.07 

438.5  nm 

3 

3.28 

+ 

0.01 

K 

10 

7469  . 

13  B  Cyg 

0.71 

0.71 

+ 

0.06* 

H,  K 

11 

0.85 

+ 

0.08* 

H,  K 

12 

7557  . 

53  a  Aql 

3.05 

3.46 

+ 

0.04 

451,  550,  800  nm 

5 

2.98 

+ 

0.14 

460.8  nm 

3 

8162  . 

5  a  Cep 

1.44 

1.52 

+ 

0.07* 

656  nm 

13 

References. — (1)  Nordgren  et  al.  1999;  (2)  Lane  et  al.  2001;  (3)  Hanbury  Brown 
et  al.  1974;  (4)  Kervella  et  al.  2003;  (5)  Mozurkewich  et  al.  2003;  (6)  Kervella  et  al. 
2004;  (7)  Shao  et  al.  1988;  (8)  Nordgren  et  al.  2001;  (9)  Boden  et  al.  2000;  (10)  Ciardi 
et  al.  2001;  (11)  Boden  &  Lane  2001;  (12)  Torres  et  al.  2002;  (13)  Vakili  et  al.  1998. 

Note.— Asterisk  (*)  denotes  measurements  we  have  corrected  to  account  for  limb- 
darkening  as  described  in  the  text. 


Fig.  1. — Comparison  between  measured  and  calculated  angular  diameters 
of  early-type  stars.  The  dashed  line  at  45°  represents  perfect  agreement,  while 
the  solid  line  shows  the  least-squares  fit  of  the  measurements  to  the  calculated 
values  with  an  additive  offset  (of  0.084  ±  0.026  mas).  The  deviations  of  (3 
Cas  and  a  Aql  are  probably  due  to  rotation. 


and  luminosity.  Interpolation  in  the  Claret  grid  (Claret  2004), 
assuming  solar  composition,  yields  a  new  mass  estimate  (e.g., 
Peterson  et  al.  2006a).  The  entire  process  typically  converges 
to  a  fraction  of  a  percent  in  mass  in  a  few  iterations. 

Once  the  low  inclination  limit  is  established  (at  co  = 
0.999  or  as  limited  by  the  ZAMS  gravity;  see  below),  we 
sample  the  visibilities  over  the  allowed  range  of  inclinations 
and  position  angles.  To  keep  calculations  as  realistic  as  possible, 
in  the  catalogs  and  the  examples  discussed  below,  the  visibil¬ 
ities  are  calculated  assuming  the  object  is  on  the  meridian,  and 
we  account  explicitly  for  its  zenith  distance  for  each  specific 
observatory. 

3.2.  The  ZAMS  Constraint 

In  the  process  of  conducting  these  simulations,  we  have 
found  that  the  deduced  polar  gravity  log  increases  mono- 
tonically  as  the  inclination  angle  decreases,  with  the  projected 
velocity  held  fixed.  Figure  2  shows  two  examples  of  this,  for 
iS  Leo  (HR  4534)  and  f  Vir  (HR  5107).  Figure  2a  illustrates 
the  case  of  /3  Leo,  whose  polar  gravity  (solid  line)  always 
exceeds  the  polar  gravity  of  the  predicted  ZAMS  (dotted  line). 
In  such  cases,  we  assume  they  are  viewed  equator-on  only  and 
are  considered  to  be  on  the  ZAMS  within  observational  error 
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Fig.  2. — Relationships  between  polar  gravity  (the  logarithm)  and  inclination 
of  (a)  (3  Leo  and  (b)  f  Vir  {solid  lines).  The  inclination  covers  the  range  from 
equator-on  to  breakup  (w  =  0.999).  The  dotted  line  shows  the  polar  gravity 
at  the  ZAMS.  /3  Leo  has  a  larger  polar  gravity  than  if  it  were  on  the  ZAMS 
with  its  predicted  mass  for  all  inclinations.  In  such  cases,  we  assume  it  is 
viewed  at  i  =  90°.  The  permitted  range  of  inclinations  for  f  Vir  is  terminated 
well  before  breakup  (co  =  0.978)  by  the  ZAMS  constraint. 

and  uncertainties  in  composition,  f  Vir,  shown  in  Figure  2b, 
illustrates  the  case  in  which  the  potential  range  of  the  incli¬ 
nations  is  terminated  before  reaching  the  inclination  corre¬ 
sponding  to  breakup,  here  taken  to  be  o)  =  0.999. 


3.3.  Statistical  Calculations 

Once  the  range  of  possible  inclinations  is  established,  is 
tabulated  over  a  grid  of  20  evenly  spaced  position  angles  and 
20  inclinations.  From  this  tabulation,  we  generate  histograms 
describing  the  probability  density  function  and  calculate  the 
maximum,  minimum,  average,  and  68.3%  (1  a)  confidence 
interval  of  V^,  found  by  splitting  the  out-of-range  probability 
equally  on  the  two  sides  (see  Yoon  et  al.  2006  for  details).  The 
assumption  here  is  that  the  rotational  axis  is  randomly  oriented 
in  space  over  the  range  of  allowed  inclinations  and  position 
angles  (e.g.,  Deutsch  1969). 

4.  RESULTS 

We  illustrate  the  range  of  results  with  four  cases  listed  in 
Table  2.  They  include  a  slow  rotator  (HR  1673),  an  intermediate 
rotator  (HR  8615),  a  moderately  fast  rotator  (HR  3974),  and  a 
fast  rotator  (HR  7740).  We  have  used  appropriate  instrument 
parameters  for  the  80  m  baseline  at  the  Navy  Prototype  Optical 
Interferometer  (NPOI;  Armstrong  et  al.  1998),  operating  at  a 
wavelength  of  500  nm.  Note  that  by  usual  standards,  these  four 
objects  would  at  first  glance  be  good  calibrator  choices  for  the 
80  m  baseline  at  the  NPOI,  with  angular  diameters  nominally 
in  the  range  of  0.5-0. 6  mas. 

To  put  the  calculations  into  context,  we  note  that  estimates 
of  inferred  angular  diameters  using  surface  brightness-an¬ 
gular  diameter  relations  are  typically  uncertain  by  a  few  per¬ 
cent.  Adopting  the  3%  uncertainty  suggested  by  Mozurkewich 
et  al.  (2003),  a  star  with  a  uniform  disk  angular  diameter  of 
9  =  0.55  mas,  measured  with  a  80  m  baseline  at  500  nm,  has 
a  visibility  =  0.63,  which  is  uncertain  by  about  ±0.02. 
The  effects  of  rotation  on  the  visibilities  can  be  substantially 
larger  than  this. 

Table  2  presents  the  results  of  the  simulations  for  these  four 
objects,  giving  V^,  and  ±  a  values  of  the  squared 

visibilities  of  each  star.  The  histograms  for  the  distribution 
of  the  y^’s  derived  from  these  calculations  are  shown  in  Fig¬ 
ure  3.  The  bin  size  for  the  histograms  is  a  tenth  of  the  maxi- 
mum-to-minimum  range  of  squared  visibility  if  the  range  is 
larger  than  2%  of  otherwise,  it  is  fixed  at  0.2%  of  the 
average  if  the  range  is  less  than  2%  of  V^. 


TABLE  2 

Examples  of  Effects  of  Rotation  on  Potential  Calibrators  (NPOI  80  m  baseline  at  X  =  500  nm) 


V  sin  (•  flBEM78‘'  p" 

HR  Name  Sp.  Type  V  B  —  V  (km  s“‘)  (mas)  (mas)  —a  +a 


1673  .  68  Eri  F2  V  5.12  0.44  10  0.606  39.99  0.6610  0.6298  0.6306  0.6304  0.6306 

3974  .  21  LMi  A7  V  4.48  0.18  148  0.546  35.78  0.6881  0.6778  0.6829  0.6783  0.6874 

7740  .  33  Cyg  A3  IV-Vn  4.30  0.11  268  0.533  21.41  0.7555  0.6054  0.6530  0.6272  0.6787 

8615  .  31  Cep  F3  III-IV  5.08  0.39  85  0.571  17.70  0.7102  0.6486  0.6576  0.6505  0.6632 


Note.— Data  represent  examples  of  the  effects  of  rotation  on  potential  calibrators  as  measured  by  the  maximum,  minimum,  weighted 
average,  and  ±  a  values  of  using  the  NPOI  (80  m  baseline  at  X  =  500  nm). 

^  Angular  diameter  0bem78  from  the  Barnes  et  al.  (1978)  relation,  for  guidance  only. 

Parallax  p  from  Hipparcos  (ESA  1997). 
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Fig.  3. — Histograms  of  squared  visibilities,  (a)  HR  1673  (68  Eri)  is  a  slow  rotator  with  v  sin  i  =  10  km  s“'.  In  spite  of  its  slow  rotation  and  narrow  range  of 
±  a,  the  histogram  displays  an  extended  tail,  albeit  of  low  probability,  as  shown  in  the  inset.  These  extreme  visibilities  result  from  the  fact  that  the  star  could  be 
seen  nearly  pole-on.  (b)  Histogram  for  HR  8615  (31  Cep),  an  intermediate-velocity  rotator  with  v  sin  i  =  85  km  s“'.  This  star  can  be  a  calibrator,  but  the  full 
range  of  squared  visibilities  should  be  properly  considered,  because  of  the  extended  tail  of  possible  C  values,  and  particularly  the  asymmetric  distribution, 
(c)  HR  7740  (33  Cyg)  is  a  fast  rotator,  with  v  sin  i  =  268  km  s"*.  As  expected,  the  rotation  significantly  affects  the  visibilities,  (d)  Moderately  rotating  star,  HR 
3974  (21  LMi),  with  v  sin  i  =  148  km  s“',  which  is  on  or  very  near  the  ZAMS  and  is  therefore  assumed  seen  equator-on.  This  can  be  a  good  calibrator,  because 
of  the  small  predicted  range  in  C,  owing  to  the  constrained  range  of  inclination. 


The  usual  assumption  that  slow  projected  rotation  assures 
no  significant  effects  from  rotation  is  not  so  reliable.  The  slow 
rotator  HR  1673  (Fig.  3a)  shows  an  extended,  low-probability 
tail  in  squared  visibility.  Even  though  the  1  a  confidence  in¬ 
terval  for  this  object  is  very  small,  there  is  a  finite  chance  of 
introducing  a  4%  error.  In  the  simulations  for  this  star,  the 
inclination  was  cut  off  at  an  inclination  of  i  =  2.45°,  with 
smaller  inclinations  pushing  the  models  below  the  ZAMS.  Even 
with  critical  rotation  excluded,  the  high  inclination  that  is  ac¬ 
cessible  can  generate  significant  uncertainties  in  the  visibilities. 

Figure  3b  shows  the  corresponding  histogram  for  HR  8615, 
which  would  normally  be  considered  a  slow  rotator  for  these 
stars,  with  a  projected  rotational  velocity  of  85  km  s  '.  The 
histograms  show  a  rather  different  story,  with  a  wing  on  the 
distribution  to  higher  squared  visibilities  amounting  to  10%  of 
the  mean  value.  Thus,  although  this  object  appears  to  be  as 
reasonable  a  candidate  for  a  calibrator  as  HR  1673,  it  should 
be  used  only  if  nothing  else  is  available. 

For  the  moderately  fast  rotator  HR  7740  (v  sin  i  =  268  km 
s“';  Fig.  3c),  we  see  the  expected  effects  on  the  predicted 


squared  visibilities.  Significant  flattening  is  predicted  at  all  in¬ 
clinations,  so  the  visibility  measurements  vary  significantly 
over  the  possible  combinations  of  inclination  and  position  an¬ 
gle.  In  spite  of  a  small  predicted  angular  diameter,  0.533  mas, 
the  range  in  visibilities  must  be  considered  if  calibrating  with 
this  star. 

Finally,  Figure  3d  shows  how  the  ZAMS  constraint  affects 
the  visibility  range  of  a  star  with  moderately  fast  projected 
rotation;  HR  8615,  with  f  sin  i  =  148  km  s  '.  This  star  appears 
to  lie  near  but  slightly  below  the  ZAMS  (which  we  interpret 
to  mean  that  it  is  on  the  ZAMS  to  observational  limits),  which 
greatly  limits  the  predicted  range  of  visibilities.  We  predict  HR 
8615  would  be  a  good  calibrator,  in  spite  of  its  significant 
projected  velocity.  The  double-horned  shape  of  the  resulting 
histogram  is  characteristic  of  stars  that  display  significant  ro¬ 
tation  and  are  near  the  ZAMS. 

From  these  examples,  it  is  clear  that  the  probability  density 
functions  for  the  squared  visibilities  are  distinctly  non-Gaus- 
sian,  and  there  is  a  significant  question  of  how  best  to  char¬ 
acterize  their  shapes.  For  the  catalogs,  we  have  decided  to  limit 
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TABLE  3 

NPOI  80  m  Baseline  (X  =  500  nm)  Catalog 


HR 

(1) 

Name 

(2) 

a  (J2000.0) 
(3) 

6  (J2000.0) 
(4) 

Sp.  Type 
(5) 

V 

(6) 

B-V 

(V) 

V  sin  i 
(8) 

P 

(9) 

(10) 

vL 

(11) 

V 

(12) 

—  o 

(13) 

+  CT 

(14) 

Q 

(15) 

E 

(16) 

82  . 

Up  And 

00  21  07.3 

-1-37  58  07 

F5  III 

5.18 

0.42 

41 

20.42 

0.7019 

0.6494 

0.6528 

0.6502 

0.6540 

A 

- 

269  . 

37/r  And 

00  56  45.2 

-1-38  29  58 

A5  V 

3.87 

0.13 

72 

23.93 

0.5870 

0.5221 

0.5297 

0.5242 

0.5330 

A 

- 

343  . 

33fl  Cas 

01  11  06.2 

-1-55  08  59 

Al  V 

4.33 

0.17 

102 

23.73 

0.6910 

0.6338 

0.6426 

0.6358 

0.6469 

A 

- 

417  . 

48cd  And 

01  27  39.4 

-1-45  24  24 

F5  IV 

4.83 

0.42 

69 

35.33 

0.6183 

0.5496 

0.5563 

0.5510 

0.5598 

A 

- 

531  . 

53x  Cet 

01  49  35.1 

-10  41  11 

F3  III 

4.67 

0.33 

61 

42.35 

0.6152 

0.6080 

0.6103 

0.6089 

0.6115 

A 

0 

Note.— Table  3  is  published  in  its  entirety  in  the  electronic  edition  of  the  PASP.  A  portion  is  shown  here  for  guidance  regarding  its  form  and  content.  Units 
of  right  ascension  are  hours,  minutes,  and  seconds,  and  units  of  declination  are  degrees,  arcminutes,  and  arcseconds.  Col.  (15):  A,  B,  and  C  characterize  the  1  a 
width  of  the  visibility  distribution  (see  text);  col.  (16):  the  characters  +,  0,  and  —  indicate  how  extended  the  tails  in  the  visibility  distributions  are  (see  text). 


the  quantitative  description  to  simply  listing  the  mean,  the  1  a 
confidence  interval,  and  the  extreme  values  of  the  squared  vis¬ 
ibilities.  Since  these  are  a  rather  minimal  set  of  parameters, 
given  the  obviously  complex  shapes  involved,  we  have  decided 
to  include  graphs  of  the  histograms  for  these  objects  along  with 
the  archived  tables. 

We  have  assumed  solar  composition  throughout  these  cal¬ 
culations.  To  see  how  sensitive  these  results  are  to  composition, 
we  have  run  the  simulations  in  a  few  representative  cases  with 
Z  =  0.03  (using  for  this  purpose  the  Padova  grid  [Girardi  et 
al.  2000],  where  Z  =  0.019  is  solar  composition).  With  the 
higher  metallicity,  the  dednced  masses  are  larger,  as  expected. 
However,  the  deduced  diameters  decrease  only  slightly,  and  in 
turn  the  visibilities  are  minimally  affected.  In  the  range  of 
diameters  enconntered  here,  the  mean  sqnared  visibilities  typ¬ 
ically  increase  by  0.3%,  with  an  occasional  case  reaching  1%. 
We  conclude  that  the  predictions  here  are  not  sensitive  to  com¬ 
position  so  long  as  it  is  not  too  far  from  solar. 

5.  CATALOGS 

The  extent  to  which  rotation  introduces  errors  in  the  cali¬ 
bration  process  depends  on  the  wavelengths  and  baselines  used. 
We  consider  three  instruments  here:  the  NPOI,  working  at 
500  nm,  the  array  at  the  Center  for  High  Angular  Resolution 
Astronomy  (CHARA;  ten  Brummelaar  et  al.  2005)  operating 
at  the  K  band,  and  the  Very  Large  Telescope  Interferometer 
(VLTI;  Glindemann  et  al.  2003),  also  at  the  K  band. 

For  the  northern  instruments,  onr  inpnt  list  of  potential  cal¬ 
ibrators  satishes  the  following  criteria:  (1)  luminosity  class  III, 
IV,  and  V  stars  in  the  color  range  B  —  V  <  0.45,  as  listed  in 
the  Bright  Star  Catalogne  (BSC;  Hoffleit  &  Warren  1991),^ 
(2)  declinations  above  —15°,  (3)  angular  size  0bem78>  estimated 
using  the  surface  brightness-angular  diameter  relation  of 
BEM78  between  0.5  and  1.0  mas,  and  (4)  distance  closer  than 
100  pc  (the  last  eliminates  significant  reddening).  We  then 
pruned  the  list  of  any  stars  with  companions  with  separations 
p  <  1"  and  magnitude  differences  Am  <  5.0,  and  of  all  con- 
hrmed  spectroscopic  binaries. 

A  separate  list  with  the  same  characteristics  but  with  a  dec¬ 


’  See  http://cdsweb.u-strasbg.fr/cgi-bin/Cat7V/50. 


lination  range  of  — 90°  <  6  <  -1-15°  was  used  for  the  VLTI 
catalog. 

The  spectroscopic  binary  (SB)  issne  is  complicated  for  the 
B,  A,  and  F  stars  and  requires  further  comment.  The  A  stars 
in  particular  are  subject  to  a  variety  of  nonradial  pulsation 
modes  (6  Scuti,  b  Delphini,  etc.)  and  inhomogeneous  surface 
composition  that,  conpled  with  rotation,  lead  to  time-variable 
structure  in  the  line  prohles.  These  have  on  numerous  occasions 
led  to  reports  of  multiple  spectra,  radial  velocity  variability, 
and  even  preliminary  SB  1  orbits  that  were  later  questioned  or 
retracted  (e.g.,  Abt  &  Levy  1974).  Unfortnnately,  once  listed, 
the  “SB”  classification  has  rarely  been  expunged  from  the  BSC. 
The  majority  of  the  A  stars  in  our  original  candidate  list  were 
characterized  as  SB.  We  have  therefore  removed  objects  re¬ 
ported  as  SB  (of  any  flavor)  only  if  a  critical  evaluation  of  the 
orbit  has  led  to  the  object’s  inclusion  in  the  Ninth  Spectroscopic 
Binary  Catalogue  (Pourbaix  et  al.  2005).®  Failing  this,  we  have 
retained  the  object  in  onr  catalog.  As  a  result,  we  expect  that 
there  will  be  binaries  lingering  in  these  lists,  and  we  repeat  the 
caution  to  always  observe  a  check  star  as  well  as  a  primary 
calibrator.  A'-band  observations  will  be  particnlarly  susceptible 
to  these  problems. 

We  have  calculated  catalogs  for  hve  specihc  configurations. 
Our  main  interest  is  in  establishing  reliable  calibrators  for  the 
NPOI,  and  three  of  the  conhgurations  are  for  the  64,  80,  and 
100  m  baselines  (AE-W7,  E6-W7,  and  E7-W7;  Armstrong  et 
al.  1998).  In  addition,  to  illustrate  the  wavelength  dependence 
of  the  results,  a  catalog  has  been  compiled  for  the  331  m 
CHARA  arm  (Sl-El;  ten  Brummelaar  et  al.  2005).  A  wave¬ 
length  of  2.22  pm  was  adopted,  typical  of  the  K  band  but  offset 
enough  from  band  center  to  avoid  Br7,  which  is  a  very  strong 
feature  in  these  objects.  Einally,  to  extend  this  to  the  Southern 
Hemisphere,  we  have  compiled  a  fifth  catalog  using  the  pa¬ 
rameters  for  the  140  m  VLTI  arm  (B3-MO;  Glindemann  et  al. 
2003),  also  at  2.22  ^tm. 

While  these  five  catalogs  are  directly  usable  for  the  instruments 
and  configurations  indicated,  they  also  provide  some  general 
guidance  as  to  the  reliability  of  individual  objects  as  calibrators. 
However,  for  instruments  and  conhgnrations  far  from  those  con- 


‘  See  http://cdsweb.u-strasbg.fr/viz-bin/Cat7V/122. 
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sidered,  we  would  be  glad  to  make  additional  runs,  given  the 
instrumental  details  (contact  the  first  author,  J.  Y). 

The  catalogs  consist  of  a  text  file  of  stars  and  their  visibility 
information,  as  shown  in  Table  2,  and  a  separate  archive  of 
postscript  files  containing  histograms:  CHARA  331  m,  NPOI 
64  m,  NPOI  80  m,  NPOI  100  m,  and  VLTI  140  m,  available  as 
electronic  files  in  the  electronic  edition  of  the  PASP.  Table  3 
shows  the  first  few  entries  of  the  NPOI  80  m  baseline  catalog, 
which  includes  the  HR  number,  name,  a  and  d  (J2000.0),  spec¬ 
tral  type,  V  magnitude,  B  —  V  color,  projected  rotational  ve¬ 
locity  (f  sin  i,  primarily  from  the  BSC),  parallax  p, 
lA,  —a,  and  +a  (in  cols.  [1]-[14],  respectively).  The  last  two 
columns,  labeled  Q  and  E,  give  a  rough  indication  of  the  quality 
of  each  object  as  a  calibrator.  In  the  column  labeled  Q,  the 


calibrator  is  classified  as  A  if  AalV^  <  0.02,  B  if  0.02  < 
A<t/\P  <  0.04,  and  C  otherwise,  where  Act  =  (+ct)  —  (— ct). 
The  last  column,  E,  classifies  a  calibrator  plus,  zero,  or  minus 
if,  respectively,  —  V^i^)IAa  <  2,  is  between  2  and  4,  or  is 
greater  than  4.  The  E  classification  is  meant  to  warn  of  cases 
in  which  one  might  be  “exposed”  to  large  error  in  V^,  even 
when  Act  is  small  and  the  probability  of  such  exposure  is  small. 
The  histograms  have  a  format  similar  to  the  figures  shown  in 
this  paper. 

This  research  is  partially  supported  by  a  grant  from  the  Naval 
Research  Laboratory  to  D.  M.  P.  We  would  like  to  thank  Dan¬ 
ielle  Kumpulanian  for  assistance  during  the  early  phases  of  this 
work. 
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